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ABSTRACT 

Discovered after the end of the Compton Gamma-Ray Observatory mission, the radio pulsar 
PSR J202H-3651 was long considered a likely counterpart of the high-energy 7-ray source 2CG 075-f-OO 
— 3EG J2021-I-3716 = GeV J2020-I-3658, but it could not be confirmed due to the lack of a contempora- 
neous radio pulsar ephemeris to fold the sparse, archival 7-ray photons. Here, we report the discovery 
of 7-ray pulsations from PSR J2021-K3651 in the 100-1500 MeV range using data from the AGILE 
satellite gathered over 8 months, folded on a densely sampled, contemporaneous radio ephemeris ob- 
tained for this purpose at the Green Bank Telescope. The 7-ray pulse consists of two sharp peaks 
separated by 0.47±0.01 cycles. The single radio pulse leads the first 7-ray peak by 0.165±0.010 cycles. 
These properties are similar to those of other 7-ray pulsars, and the phase relationship of the peaks 
can be interpreted in the context of the outer-gap accelerator model for 7-ray emission. Pulse-phase 
resolved images show that there is only one dominant source, AGL J2020.5-f 3653 — PSR J2G21-I-3651, 
in the region previously containing confused sources 3EG J2021-)-3716 and 3EG J2016+3657. 
Subject headings: gamma-rays: observations — pulsars: individual (PSR J2G21+3651) 



1. INTRODUCTION 

Rotation-powered pulsars are the dominant class 
of high-energy (> 100 MeV) 7-ray sources in the 
Galaxy. The Energetic Gamma-Ray Experiment Tele- 
scope (EGRET) on the Gompton Gamma-Ray Observa- 

■ tory (CGRO), which operated between 1991 and 2000, 
detected si x pulsars, and t hree others with less confi- 

■ dence. (See lThompsonll2004l for a review of these EGRET 
results.) It is widely assumed that a large number of 
unidentified 7-ray sources in the Galactic plane are also 
pulsars. Some of these are positionally coincident with 
radio pulsars, but do not have enough photons to recover 
the pulse period, while others may be radio quiet. An 
ephemeris from radio or X-ray observations is generally 

' needed to fold the sparse 7-ray source photons over a 
span of months to years in order to build up a statisti- 
cally significant pulsed light curve. In several cases, a 
radio pulsar candidate was only discovered after the end 
of the CGRO mission, so that the needed contemporane- 

' ous ephemeris does not exist. The typical "timing noise" 
of young pulsars causes a phase ephemeris to fail when 
extrapolated over months to years. 

Some of the still-unidentified Gal actic 7-ray sources 
date from the 1975 GOS B mission (|Swanenburg et all 
119811 ). The target of this work is the GOS B source 
2CG 075-j-OO located in Cygnus. EGRET possibly re- 

^ Columbia Astrophysics Laboratory, Columbia University, New 
York, NY 10027. 

2 INAF/IASF-Milano, 1-20133 Milano, Italy 

^ Department of Physies, West Virginia University, Morgantown, 
WV 26501. 

^ Department of Physics and Astronomy, Barnard College, 
Columbia University, New York, NY 10027. 

^ National Radio Astronomy Observatory, Charlottesville, VA 
22903. 

" Eureka Scientific, Inc., Oakland, CA 94602. 

INAF/IASF-Roma, 1-00133 Roma, Italy 
* Dipartimento di Fisica, Universita "Tor Vergata" , 1-00133 
Roma, Italy. 



solved 2CG 075-t-OO into two sources, 3EG_J2021+3716 
and the weaker 3EG J2016-I-3657 (|Hartman et al.l 
19991). Multiw a yeleng th st udy of these region s by 
Mukherjee et all (|200d l and iHalpern et all (|2001h un- 
covered possible counterparts for 3EG J2016-I-3657, 
but no obvious candidate for 3EG J2021-F3716. A 
smalle r error region was derived by iLamb fc MacombI 
(1997*) by restricting analysis to energies above 1 GeV, 
yielding a source GeV J2020H-3658 that overlapped 
with 3EG J2021-h3716. An ASGA observation of 
GeV J2020-I-3658 fo und two compact X-ray sources 
([Roberts et al.l [200lh . one of which was discovered to 
cont ain a new, energetic radio pulsar, PSR J2021-f3651 
(.Roberts et al.|[2002,) . Its spin period P = 0.103 s, char- 
acteristic age Tc = P/2P — 17 kyr, and spin-down lumi- 
nosity E — 3.4 X 10^^ ergs s~^ are in the range of other 
young pulsars that were detected by EGRET, but its ini- 
tial distance estimate of ~ 12 kpc implied an unusually 
large 7-ray efficiency (see §3.3) if PSR J2021-I-3651 is re- 
sponsible for 3EG J2021-I-3716. It was not possible to ex- 
trapolate the radio ephemeris backward to search for the 
pulsar in EGRET 7-ray photons. By searching in a range 
around the expected P and P, iMcLaughlin fc Corded 
(|2004f l possibly detected pulsations, but only in two out 
of eight EGRET viewing periods. Thus, the proposed 
identification of PSR J202H-3651 with 3EG J2G21+3716 
and GeV J2020-I-3658 could not be confirmed, until now. 

The AGILE (Astro-rivelatore Gamma a Immagini 
LEggero) mission of the Italian Space Agency (ASI) 
was launched on 2 007 April 23 into a low-Earth orbit 
(jTavani et al.ll2008f ). AGILE will allow a number of new 
7-ray pulsars to be found by folding on contemporane- 
ous pulsar ephemerides that are being constructed by 
our group and others using several radio telescopes. For 
the analysis of PSR J2021-I-3651, we have been using 
the NRAO Green Bank Telescope (GBT). In this Let- 
ter, we apply the ephemeris of PSR J2021+3651 to data 
obtained by AGILE during the first year of the mission. 
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2. OBSERVATIONS AND ANALYSIS 

The AGILE GRID (Gamma-Ray I maginR Detector) 
comprises a Silicon- Tu ngsten tracker ("Barbielli ni et alj 
120011: iPrest et al.ll2003D and a CsI(Tl) mini-calorimeter 
( Labanti et al.l |2006[ ). These provide effective area of 
~ 500 cm^ in the 30 MeV-50 GeV range over a ^ 2.5 sr 
field of view. We used GRID observations of the Cygnus 
region during five periods between 2007 November 2 and 
2008 June 30 that sum to 103 days of observation when 
PSR J2021-h3651 was within 40° of the pointing direc- 
tion. GRID photons are time-tagged with a precision of 
^ 1 /is, and the arrival times were transformed to the 
solar system barycentcr in Barycentric Dynamical Time 
using GPS positioning of the spacecraft and the Chan- 
dra measured position of the pulsar, (J2000.0) R.A. = 
20'^21^°05.46^ decl. = -h36°51'04.8", uncertainty - 0.5" 
(|Hessels et al.l [20041. The accuracy of this procedure is 
~ 10 /iS, and the absolute times were verified to a pre- 
cision of ^ 400 pLS using observat ions of the Crab pulsar 
and others (jPellizzoni et al.ll2008D . 

A phase-connected rotational ephemeris for 
PSR J2021+3651 contemporaneous with all of the 
AGILE observing periods was derived from 36 observa- 
tions obtained at the GBT between 2007 September 28 
and 2008 June 30. The pulsar was obs erved with the 
GBT Pulsar Spigot (jKaplan et al.ll2005D . yielding total 
power samples every 81.92 iis in each of 768 frequency 
channels spanning a bandwidth of 600 MHz centered 
on 1950 MHz. Each observation usually lasted for 5 
minutes, from which we derived a time of arrival with 
typical uncertainty of 0.2 ms. We used the TEMPO^ 
timing software to obtain the ephemeris using standard 
methods. Like many young pulsars, PSR J2021-f3651 
exhibits significant rotational instability. We parame- 
terized this by fitting for the rotational frequency and 
its first six derivatives. This timing solution had small 
unmodeled residual features {xl = 2.5), with a post-fit 
rms of 0.33 ms and maximum of 0.9 ms. This precision 
is comparable to the absolute accuracy of the AGILE 
timing, and ensures that the effective time resolution of 
the folded 7-ray light curve is ^ 0.01 P. The observed 
radio dispersion measure (DM = 369.3 ± 0.3 pc cm~^, 
determined with the addition of a few observations 
centered at 1550 MHz) was used to correct the times 
of arrival of the radio pulse to infinite frequency (with 
resulting la uncertainty of 0.6 ms) for absolute phase 
comparison with the 7-ray pulse. 

In choosing the sky region around the pulsar posi- 
tion for extraction of 7-ray photons, consideration of the 
point-spread function (PSF) of the instrument, which in- 
creases with decreasing energy, must be weighed against 
the high diffuse 7-ray background in this region, which 
also increases with decreasing energy, as well as possi- 
ble contamination from nearby point sources. We find 
that the pulsed power is maximized if we use an extrac- 
tion region that is comparable to the 68% containment 
radius, which is 3.5 ° at 100 MeV and 1.2° at 400 MeV 
(jTavani et al.ll2008f ). We used r = 1.2° for E > 400 MeV, 
and r = 1.2° {400/ E)^/'^ for E < 400 MeV. In addition, 
we use only the "G" class events, which are identified 
with good confidence as photons, and exclude the "L" 

^ See |http: / / www.atnf.csiro.au/research/pulsar/tempol 
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Fig. 1. — Folded 7-ray light curves (histograms) of 
PSR J2021+3651 in two energy bands, and their sum. Two cycles 
are plotted, with 40 bins per cycle. The solid line in the bottom 
panel is the mean radio pulse at 1950 MHz from the GBT. Phase 
lags are P2-P1=0.47 ± 0.01 cycles and Pl-R=0.165 ± 0.010 cycles. 

class which suffer much higher particle background con- 
tamination. This selection limits the effective area to 
~ 250 cm^ at 100 MeV. 

3. RESULTS AND INTERPRETATION 

3.1. Pulsed Light Curve 

Figure [T] shows the resulting folded light curve of 
PSR J2021+3651 in two energy bands, 100 - 300 MeV 
and 300 — 1500 MeV. The pulse shape consists of two 
sharp peaks, labeled PI and P2, with separation of 
0.47 ± 0.01 cycles (consistent with the separation in 
the tentative detections rep orted in EGRET data by 
iMcLaughlin fc Corded [20041 ). This structure is similar 
to that seen with much better signal-to-noise ratio in 
other pulsars. For comparison, peak separations are 0.38, 
0.43, and 0.49 for the Crab, Vela, and Geminga, respec- 
tively. As in the other 7-ray pulsars, the radio pulse 
of PSR J2021-I-3651 leads the 7-ray PI by a significant 
phase, in this case 0.165±0.010 cycles. In Vela, the radio 
leads by 0.13 cycles. The Crab pulsar has multiple broad 
radio peaks that continue to high frequency, and coincide 
with higher-energy emission from infrared through X-ray. 
However, if we identify the narrow low-frequency compo- 
nent (LEG: Moffett & Hankins 1996, 1999) as the usual 
polar-cap core emission analogous to other pulsars, the 
7-ray~radio lag is 0.11 cycles. Alternatively, the Crab's 
even lower frequency "precursor" pulse, which leads the 
7-ray PI by 0.056 cycles, may be the core component. 
Even in the case of PSR B1706-44, where AGILE re- 
solved two peaks separated by 0.24 cycles from the dom- 
inant " bridge" emission, the radio pulse leads PI by 0.21 
cycles ()Pellizzoni et al.ll2008l ). 
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Fig . 2. — Reproduction of Figure 3 of IRomani fc Yadigaroglul 
comparing 7-ray peak separation A with 7-ray-radio lag 5 
for CGRO pulsars (filled circles). The 5 of the Crab pulsar (0531) 
is uncertain due to the ambiguous identification of its core radio 
pulse. The location of PSR J2021+3651 is marked wit h an asterisk, 
as is the closely situated PSR B1951+32 ( Ra manamurthv et al.l 
119951) . Model predictions (solid lines) for different magnetic axis 
inclination angles a are similar, while viewing angle ( increases 
along the curves from lower right to upper left. 

The outer-gap mod e l of 7-ray pulse shapes of 
IRomani fc Yadigaroglul (|1995D predicts an inverse de- 
pendence of the 7-ray peak separation on the 7-ray- 
radio lag, with the principal controlling parameter be- 
ing viewing angle C h'om the rotation axis. Large C re- 
sults in large peak separation, and small 7-ray~-radio lag. 
PSR J2021-)-3651 is roughly in accord with this predic- 
tion (Fig. [2]). In this picture, pulsars that have peak 
separations close to 0.5, such as Vela, Geminga, and 
PSR J2021+3651, have C close to 90°. Such an angle 
for PSR J2021-I-3651 is also supported by the modcl- 
ing of its X-ray pu l sar wind nebula (PW N) G75. 2-1-0.1 
(|Hessels et al.ll2004l: lymT Etten et al.ll2008| l. which can be 
understood as a nearly edge-on torus, indicating a pul- 
sar rotation axis nearly perpendicular (86°) to the line 
of sight. In contrast, PSR B1706— 44 has narrow 7-ray 
peak separati on of 0.24 cycles, and ( = 53° from its X-ray 
torus model. iNg fc Romanil (|2008f ) summarize the X-ray 
measurements and their implications for 7-ray pulse pro- 
files. 

3.2. Position, Flux, and Spectrum 

In order to separate the contribution of 
PSR J2021-I-3651 from other possible 7-ray sources 
within its PSF, we made images corresponding to 
on-pulse and off-pulse intervals of its light curve (Fig.[3|). 
The on-pulse image includes six phase bins on each 
of PI and P2, while the off-pulse image comprises the 
remaining 28 bins. A difference image was created 
by subtracting the off-pulse image from the on-pulse 
image, scaled by the number of phase bins. As expected, 
differencing removes all diffuse background and sources 
except for PSR J2021-I-3651. We measure the position 
of the resulting source AGL J2020.5-I-3653 as (J2000.0) 
R.A. = 20h20'"30^ decl. = -h36°53.7', which is only 
0.12° from the X-ray position, well within the current 
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Fig. 3. — Images of the AGILE field centered on 
PSR J2021-I-3651 (cross) in the 100-1500 MeV range. Pulse- 
phase resolved images were created as described in §3.2, and 
PSR J2021+3651 was isolated via difference imaging. North is 
up and east is to the left. The diameter of the field displayed is 
20°, and the pixel size is 0.125°. Gaussian smoothing of cr = 1.5 
pixels has been applied. Each image is linearly scaled to its highest 
point (white), while the inset stretches the scale to show the level 
of fluctuations in the off-pulse image in the vicinity of the pulsar. 



AGIL E systematic uncertainty of 0.25° (jTavani et al.l 
[2OOI . We performed a standard spectral analysis on 
the total source photons including exposure correction 
and modelling of the diffuse Galactic background. The 
energy range 100—1000 MeV is fitted with a power 
law of photon index P = 1.86 ± 0.18 and a flux of 
(46 ± 6) X 10^* photons cm^^ s^^, while a turndown is 
seen above 1.5 GeV (Fig. SJ. 

Using the counting statistics of the off-pulse image, we 
estimate the level of fluctuation that would correspond 
to a 3cr point source at the location of PSR J2021-I-3651, 
noting that there is no source present at this level. Scal- 
ing from the pulsar flux, this corresponds to upper limits 
of 1.0 X 10^^ photons cm~^ for the off-pulse flux from 
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Fig. 4.— Photon spectrum of AGL J2020. 5-1-3653 with power-law 
fit of F = 1.86 to the first three points. 
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PSR J2021+3651, and 1.7 x 10"'^ photons cm"^ g-i ^j. 
a steady point source within 1.2° of PSR J2021+3651, 
such as the PWN G74.9+1.2 = CTB 87 (of unknown 
spin parameters) or the blazar B2013+ 370, both located 
in the error circle o f 3EG J2016+3657 (jMukheriee et all 
l2000HHalpern et al... 2001 1. Since blazars are episodic 7- 
ray emitters, it cannot be ruled out that B2013+370 was 
more active during the EGRET era. However, we are 
not confident that 3EG J2016+3657 was a source distinct 
from 3EG J2021+3716. T he former is not included in th e 
revised EGRET catalog of lCasandiian fc Grenied ||2008D . 
which instead lists only one source, EGR J2019+3722. 
The positions of 3EG J2021+3716 and EGR J2019+3722 
are both inconsistent with PSR J2021+3651 at their 99% 
confidence contours, while AGILE sees one dominant 
source that is clearly attributed to PSR J2021+3651 by 
phase-resolved imaging. 

3.3. Gamma-Ray Efficiency 

Previous authors have noted that the 7-ray lumi- 
nosity of 3EG J2021-I-3716 is a large fraction 77^ of 
the spin-down luminosity of PSR J2021-t-3651 if placed 
at the nomina l DM distance of 12 kpc using the 
ICordes fc Lazid (l2002lj electron-d ensity model of the 
Galaxy. [Roberts et all (|2002f ) and iHessels et all (|2004l ) 
calculated, assuming l/47r beaming, that r]^ = F(f/E « 
0.18 d?2: where F is the flux in the 100 MeV-10 GeV 
range. Such high efficiency deviates from the trend of 
other 7-ray pulsars, which follow an approximate rj-y oc 
^-1/2 relation (jThompsonll2004[ ). Pulsars of similar E 
to PSR J2021+3651 have r]^ ~ lO'^ - IQ-^. This con- 
trast casts some suspicion on the accuracy of the DM dis- 
tance, but in addition, consideration of the X-ray proper- 
ties of PSR J2021-I-3651, including the thermal spectrum 
of the neutron s t ar an d the size of its PWN torus, led 
Ivan Etten et al.l ()2008[ ) to favor a distance as small as 
3-4 kpc. 

Our results do not bear directly on the dis- 
tance question, while the inferred 7-ray efficiency of 
PSR J2021-I-3651 remains exceptional. The AGILE flux 
above 100 MeV is « 2.4 x 10~^° ergs cm~^ s~^ yield- 
ing Tj^ « 0.10(^221 somewhat smaller than previous es- 
timates b ecause of the cutoff ab o ve 1. 5 GeV. In the 
models of iRomani fc Yadigaroglul ()1995l ). the beaming 



of PSR J2021-t-3651 may be more severe than in other 
7-ray pulsars if its nearly edge-on X-ray torus implies a 
large magnetic inclination angle a. 

The TeV source MGRO J2019-h37 (|Abdo et all 
I2007allb ) of diameter 1.1° ± 0.5° centered at (£,6) = 
(75.0°, 0.2°) may be an extension of G75. 2-1-0.1, simi- 
lar to many other TeV PWNe that are spatially off- 
set from their smaller X-ray counterparts. The flux of 
MGRO J2019-h37, ^ 5.6 x lO'^^ ergs cm^^ g-i, repre- 
sents « 0.03 d?9 E of PSR J2021-t-3651, similar to other 
TeV PWNe (|Gallant et al.l[2008il . 

4. CONCLUSIONS 

After the demise of GGRO, several candidate pulsar 
identifications for EGRET sources could not be con- 
firmed until it became possible again to obtain contempo- 
raneous 7-ray and radio timing. Using the extensive ex- 
posure of A GILE on the Cygnus region, we detect 7-rays 
in the 100-1500 MeV range from the first post-EGRET 
pulsar, PSR J2021-|-3651=AGL J2020.5-I-3653. Its pulse 
shape is similar to other 7-ray pulsars, with two sharp 
peaks separated by slightly less than 0.5 cycles. This 
morphology, as well as the lag between the radio and 
first 7-ray peak, are in accord with simulations of outer- 
gap accelerators viewed at large angles from the pulsar 
rotation axis, which may help to explain its high 7-ray ef- 
ficiency. PSR J2021-I-3651 is the only 7-ray point source 
detected in the region of 2CG 075-1-00. Timing noise of 
young pulsars such as PSR J202H-3651 and other candi- 
date 7-ray pulsars demands frequent radio observations 
to be able to discover and resolve the finest structure in 
7-ray light curves that will be accumulated over months 
to years. Independent discovery of pulsed 7-rays from 
PSR J2021-I-3651 by the Fermi Gamma-ray Space Tele- 
scope will also be reported shortly. 
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